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BLIND SOURCE SEPARATION OF PULSE OXIMETRY SIGNALS 

BACKGROUND OF THE INVENTION 
[01 ] The present invention relates to the processing of signals obtained from 
a medical diagnostic apparatus such as a pulse oximeter using a blind source separation 
technique to separate the obtained data without prior knowledge of its magnitude or 
frequency into data corresponding to the desired physiological data and the undesired 

interference sources. 

[02] A typical pulse oximeter measures two physiological parameters, 
percent oxygen saturation of arterial blood hemoglobin (Sp0 2 or sat) and pulse rate. Oxygen 
saturation can be estimated using various techniques. In one common technique, the 
photocurrent generated by the photo-detector is conditioned and processed to determine the 
ratio of modulation ratios (ratio of ratios) of the red to infrared signals. This modulation ratio 
has been observed to correlate well to arterial oxygen saturation. The pulse oximeters and 
sensors are empirically calibrated by measuring the modulation ratio over a range of in vivo 
measured arterial oxygen saturations (SaO z ) on a set of patients, healthy volunteers, or 
animals. The observed correlation is used in an inverse manner to estimate blood oxygen 
saturation (SpO z ) based on the measured value of modulation ratios of a patient. The 
estimation of oxygen saturation using modulation ratios is described in U.S. Patent No. 
5,853,364, entitled "METHOD AND APPARATUS FOR ESTIMATING 
PHYSIOLOGICAL PARAMETERS USING MODEL-BASED ADAPTIVE FILTERING", 
issued December 29, 1998, and U.S. Patent No. 4,91 1,167, entitled "METHOD AND 
APPARATUS FOR DETECTING OPTICAL PULSES", issued March 27, 1990. The 
relationship between oxygen saturation and modulation ratio is further described in U.S. 
Patent No. 5,645,059, entitled "MEDICAL SENSOR WITH MODULATED ENCODING 
SCHEME," issued July 8, 1997. Most pulse oximeters extract the plethysmographic signal 
having first determined saturation or pulse rate, both of which are susceptible to interference. 

[03] A challenge in pulse oximetry is in analyzing the data to obtain a 
reliable measure of a physiologic parameter in the presence of large interference sources. 
Prior art solutions to this challenge have included methods that assess the quality of the 
measured data and determine to display the measured value when it is deemed reliable based 
upon a signal quality. Another approach involves a heuristic-based signal extraction 



technology, where the obtained signals are processed based on a series of guesses of the ratio, 
and which require the algorithm to start with a guess of the ratio, which is an unknown. Both 
the signal-quality determining and the heuristic signal extraction technologies are attempts at 
separating out a reliable signal from an unreliable one, one method being a phenomenological 
5 one and the other being a heuristic one. 

[04] On the other hand, a problem encountered in such disciplines as 
statistics, data analysis, signal processing, and neural network research, is finding a suitable 
representation of multivariate data. One such suite of methods is generally known as 
Independent Component Analysis (ICA), which is an approach to the problem of Blind 

10 Source Separation (BSS). 

[05] In general terms, the goal of blind source separation in signal 
processing is to recover independent source signals after they are linearly mixed by an 
unknown medium, and recorded or measured at N sensors. The blind source separation has 
been studied by researchers in speech processing or voice processing; antenna array 
1 5 processing; neural network and statistical signal processing communities (e.g. P. Comon, 
"Independent Component Analysis, a New Concept?", Signal Processing, vol. 36. no. 3, 
(April 1994), pp. 287-314, "Comon") and applied with relative degrees of success to 
electroencephalogram data and functional MRI imaging. 

[06] Comon defined the concept of independent component analysis as 
20 maximizing the degree of statistical independence among outputs using "contrast" functions 
of higher-order cumulants. Higher-order statistics refer to the expectations of products of 
three or more signals (e.g. 3 rd -order or 4 th -order moments), and cumulants are functions of 
the moments which are useful in relating the statistics to those of the Gaussian distribution. 
The 3 rd -order cumulant of a distribution is called a skew, and the 4 th -order cumulant is the 
25 kurtosis. A contrast function is any non-linear function which is invariant to permutation and 
scaling matrices, and attains its minimum value in correspondence of the mutual 
independence among the output components. In contrast with decollation techniques such 
as Principal Component Analysis (PCA), which ensures that output pairs are uncorrelated, 
ICA imposes the much stronger criterion that the multivariate probability density function of 
30 output variables factorizes. Finding such a factorization requires that the mutual information 
between all variable pairs go to zero. Mutual information depends on all higher-order 
statistics of the output variables while decorrelation normally only takes account of 2nd-order 
statistics. 
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107] While the general use of IC A as a means of blindly separating 
independent signal sources is known, the method poses unique challenges to its 
implementation in pulse oximetry. For instance, the mixture signals may not be exactly a 
linear combination of the pulse signal and sources of interference. Also, most ICA 
5 techniques are based on fourth-order cumulants, as the signals and noise commonly 

encountered in communications have zero third-order cumulant (skew), and cumulants of 
higher than fourth order are difficult to estimate accurately. 

[08] Several ICA methods are known for separating unknown source 
signals from sets of mixture signals, where the mixture signals are a linear combination of the 
10 source signals. As used in pulse oximetry, the mixture signals refer to signals measured at 
multiple wavelengths. Source components refer to the desired physiologic data including 
signals corresponding to the plethysmography signal obtained at multiple wavelengths in 
addition to undesired interference data, which may be caused by motion, light interference, 
respiratory artifacts, and other known sources of errors in pulse oximetry. 
15 l09 ] There is therefore a need to apply blind source separation techniques to 

the field of pulse oximetry to be able to deterministically separate a source signal from 
various interference sources. 

BRIEF SUMMARY OF THE INVENTION 
[10] The present invention is directed towards a method and apparatus for 
20 the application of Blind Source Separation (BSS), specifically Independent Component 

Analysis (ICA) to pulse oximetry. ICA refers to any one of several methods for separating 
unknown source signals from a set of "mixture" signals, which are linear combinations of the 
source signals. These methods may use estimates of the second- and higher-order joint 
statistics of the mixture signals and separate the sources by seeking to minimize the mutual 
25 information of the outputs of separation. In pulse oximetry, the signals measured at different 
wavelengths represent the mixture signals, while the plethysmography signal, motion 
artifact, respiratory artifact and instrumental noise represent the source components. 

[11] In one embodiment the BSS is carried out by a two-step method 
including PCA and a higher-order decollation. In the first step, the method uses PCA as a 
30 preprocessing step, and in a second step, the principal components are then used to derive the 
independent components and the desired physiological parameters. The PCA is performed to 
transform the data to have zero second-order correlation before higher-order decorrelation. 
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[12] In one aspect of the method of the present invention, data 
corresponding to a plurality of signals measured at a plurality of wavelengths are first 
obtained. Next, the data are processed to obtain a plurality of principal components, where in 
one embodiment the principal components are obtained by decorrelating the data (to 
5 minimize the cross-correlation between the signals from different wavelengths), and 

normalizing the decollated data. Next, the principal components are processed to obtain a 
plurality of independent components, wherein a matrix of the plurality of signals corresponds 
with a matrix product of a matrix of the plurality of independent components and a matrix of 
mixing coefficients. In one embodiment, the independent components are obtained by 
1 0 higher-order decorrelation of the principal components, and where the higher-order 

decorrelation of the principal components is achieved by minimizing a function of the higher- 
order cross-correlation of the data or equivalent^ by maximizing a function of the higher- 
order cumulants of the plurality of mixture signals. Since the skew of the time-derivative of 
the pulse signal is generally much greater in magnitude than that of interference, performance 
15 of the ICA may be enhanced by using a "contrast" function that was derived from the third- 
order cumulants of the derivatives of the signals. 

[13] In an aspect of the method of the present invention directed towards a 
pulse oximeter measuring signals at multiple wavelengths, a first independent component 
corresponds with a plethysmography signal, a second independent component corresponds 
20 with the interference sources, and sat may be determined from a ratio of mixing coefficients 
from the mixing matrix. In pulse oximetry, the technique provides the advantage of 
extracting the plethysmographic signal in the presence of large motion interference and 
especially without requiring prior knowledge of saturation or pulse rate. 

[14] For a further understanding of the nature and advantages of the 
25 invention, reference should be made to the following description taken in conjunction with 
the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[15] Fig. 1 is a block diagram of an exemplary pulse oximeter. 
[16] Fig. 2 is a flow chart of an embodiment of the method of the present 

30 invention. 

[17] Fig. 3 is a graph showing a typical pulse oximetry signals at two 

wavelengths. 
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[18] Fig- 4 is a graph showing the typical pulse oximetry signals at two 
wavelengths after PCA processing. 

[19] Fig- 5 is a graph showing a typical pulse oximetry signals at two 

wavelengths after IC A processing. 
5 [201 Fig. 6 is a graph of signals of Fig. 3 plotted against one another. 

[211 Fig. 7 is a graph of the principal components of Fig. 3 plotted against 

one another. 

[22J Fig. 8 is a graph of the independent components of Fig. 3 plotted 
against one another. 

I o DETAILED DESCRIPTION OF THE INVENTION 

[23] The present invention is directed towards the application of Blind 
Source Separation (BSS), specifically Independent Component Analysis (ICA) to pulse 
oximetry. ICA refers to any one of several methods for separating unknown source signals 
from a set of "mixture" signals, which are linear combinations of the source signals. The 

1 5 ICA method as embodied by the present invention, uses estimates of the second- and higher- 
order joint statistics of the mixture signals and separates the sources by seeking to minimize 
the mutual information of the outputs of separation. In pulse oximetry, the signals measured 
at different wavelengths represent the mixture signals, while the plethysmographic signal, 
motion artifact, respiratory artifact and instrumental and environmental noise represent the 

20 source components. 

[24] In one embodiment, the BSS is carried out by a two-step method 
including an ICA. In the first step, the method uses Principal Component Analysis (PCA) as 
a preprocessing step, and the Principal Components are then used to derive sat and 
Independent Components, where the Independent Components are determined in the second 

25 step. Before describing the BSS methods of the present invention an example of a pulse 
oximeter, which may be configured to practice the method of the present invention is 
described below. 

[25] FIG. 1 is a block diagram of one embodiment of a pulse oximeter 
implementing the present invention. Light from light source 1 10 passes into patient tissue 
30 112, and is scattered and detected by photodetector 1 14. A sensor 100 containing the light 
source and photodetector may also contain an encoder 1 16 which provides signals indicative 
of the wavelength of light source 1 10 to allow the oximeter to select appropriate calibration 
coefficients for calculating oxygen saturation. Encoder 1 16 may, for instance, be a resistor. 
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[26] Sensor 1 00 is connected to a pulse oximeter 1 20. The oximeter 
includes a microprocessor 122 connected to an internal bus 124. Also connected to the bus is 
a RAM memory 126 and a display 128. A time processing unit (TPU) 130 provides timing 
control signals to light drive circuitry 132 which controls when light source 1 10 is 
5 illuminated, and if multiple light sources are used, the multiplexed timing for the different 
light sources. TPU 130 also controls the gating-in of signals from photodetector 1 14 through 
an amplifier 133 and a switching circuit 134. These signals are sampled at the proper time, 
depending upon which of multiple light sources is illuminated, if multiple light sources are 
used. The received signal is passed through an amplifier 136, a low pass filter 138, and an 
10 analog-to-digital converter 140. The digital data is then stored in a queued serial module 
(QSM) 142, for later downloading to RAM 126 as QSM 142 fills up. In one embodiment, 
there may be multiple parallel paths of separate amplifier filter and A/D converters for 
multiple light wavelengths or spectrums received. 

[27] Based on the value of the received signals corresponding to the light 
15 received by photodetector 1 14, microprocessor 122 will calculate the oxygen saturation using 
various algorithms. These algorithms require coefficients, which may be empirically 
determined, corresponding to, for example, the wavelengths of light used. These are stored in 
a ROM 146. The particular set of coefficients chosen for any pair of wavelength spectrums is 
determined by the value indicated by encoder 1 16 corresponding to a particular light source 
20 in a particular sensor 100. In one embodiment, multiple resistor values may be assigned to 
select different sets of coefficients. In another embodiment, the same resistors are used to 
select from among the coefficients appropriate for an infrared source paired with either a near 
red source or far red source. The selection between whether the near red or far red set will be 
chosen can be selected with a control input from control inputs 154. Control inputs 154 may 
25 be, for instance, a switch on the pulse oximeter, a keyboard, or a port providing instructions 
from a remote host computer. Furthermore any number of methods or algorithms may be 
used to determine a patient's pulse rate, oxygen saturation or any other desired physiological 
parameter. One such method, namely Blind Source Separation, is described below. 

[28] Blind Source Separation refers to the separation of signals given only 
30 linear combinations of those signals, such that: 
x(t) = As(t) 

where x(t) is a matrix of a set of observed signals (mixed signals), xi(t) 

....x n (t), 

A is an unknown mixing matrix, 
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and s(t) is a set of source signals Si(t) . ..s m (t), 
assumed to be statistically independent, i.e. 

m 

p(s) = p(si,..,s m )= J~[/w(.sO 
1=1 

where p(s) is the probability distribution function of s. 
5 [29] As described above, in pulse oximetry, the mixture signals correspond 

with signals obtained by a pulse oximeter sensor, which include both the desired signal and 
the undesired noise components. In one embodiment of the method of the present invention, 
the mixture signals are first preprocessed using PCA to transform the mixture signal to 
principal components. To more fully separate the signal and the noise, the data are further 
10 processed: in mathematical terms, the data are rotated. In other words, the ICA processing 
includes a combination of PCA and rotation. 

[30] A criterion for determining the degree of signal-noise separation is 
statistical independence, as described above. However, since the probability distributions are 
not known, the challenge of an ICA algorithm becomes the measurement of statistical 
1 5 independence. A measure of statistical independence is the degree of mutual information, 
such that by minimizing the degree of mutual information between sets of data, independent 
components can be determined. Algorithms determining the mutual information are 
generally too complicated for a direct solution of its (i.e. mutual information) minimum, and 
thus they lend themselves best to iterative methods. For example, one possible approach 
20 would be to search for coefficients of the mixing matrix A that would lead to statistical 
independence (by minimizing the data set's mutual information). One could heuristically 
sweep through a large range of angles about which to rotate the principal components, which 
would yield an independent set of data, but this approach would be excessively time- 
consuming. 

25 [31 ] Thus the inventor of the present invention proposes separating the data 

by performing higher-order decorrelation of the data, or by removing the higher-order 
correlation of the data obtained from the mixture signals. Thus, the BSS-based method of the 
present invention: (1) uses PCA to find uncorrelated components, and (2) separates the data 
by removing higher-order correlation to find the independent components corresponding with 

30 the desired signal source(s) and the undesired noise sources. And as used herein, higher- 
order correlations are higher than second-order correlations, such as third-order and fourth- 
order cross-correlations. In one embodiment, independence is approximated by minimizing 
the sum of the squares of the third-order correlations r xxy and r xyy , e.g., 
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r m = 2>| l * and r *rr «IW [where x and v have zero mean] 

[32] Alternately, independence is approximated by minimizing the sum of 
the squares of the fourth-order cross-cumulants. 

[33] In certain embodiments of the present invention, separation is achieved 
by maximizing the sum of the squares of the higher-order cumulants, which is equivalent to 
minimizing the sum of the squares of the higher-order cross-cumulants. In these and other 
embodiments, the second-order decorrelation and higher-order decorrelation may be effected 
simultaneously through an iterative, adaptive process. 

[34] The advantage of achieving separation by higher-order decorrelation is 
that it enables direct formulas and simple algorithms for the separation of data into its 

independent components. 

[35] Fig. 2 is a flow chart 200 of an embodiment the method of the present 
invention as applied to signals obtained by a pulse oximeter sensor. First, (step 210) a 
plurality of signals measured at various wavelengths are obtained by a pulse oximeter sensor. 
In a typical pulse oximeter, emitting optical energy at two wavelengths, the photocurrent 
generated by the photo-detector is conditioned and processed to determine the modulation 
ratios of the red and infrared signals. The example of a two-wavelength pulse oximeter is for 
illustration purposes only, and is not meant to limit the scope of the present invention. The 
detected photocurrents include the mixture signals, where the mixture signals include the 
plethysmographic signal, motion artifact, respiratory artifact and instrument noise. An 
example plot of the measured photocurrent is shown on Fig. 3. Fig. 3 shows a plot 300 of the 
photocurrent vs. time for measurements obtained at 595 ran (310) and at 665 nm (320). As 
can be seen from this figure (Fig. 3), both signals (310 and 320) show a low amplitude 
section and a high amplitude section, where the high amplitude section corresponds to the 
signals measured while the sensor is moving. 

[36] Next, (step 220) the mixture signals are processed using PCA analysis 
to obtain two principal components. Here, the PCA processing results in the determination of 
two principal components, since the original mixture signals were obtained at two different 
wavelengths. Embodiments of the present invention are not limited to two mixture signals or 
two principal components. The embodiments of the present invention are directed to the 
decomposition of a matrix of a set of observed signals (mixed signals) into a set of source 
signals and a mixing matrix, as set forth above. However, describing more than two 
dimensions and visualizing more dimensions is at best difficult to visualize, hence the 
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description provided herein is kept to a two-dimensional one. In one embodiment, the 
Singular Value Decomposition algorithm is used to obtain the principal components. In one 
alternate embodiment, the data are multiplied by the inverse of the square root of the 
covariance matrix. In another embodiment directed to a two wavelengths approach, for each 

5 possible pair of wavelengths, the data are rotated by the angle of the best linear fit between 
those two signals. Fig. 4 shows a plot 400 of the two principal components. As can be seen 
from this figure (Fig. 4), the first principal component 410 corresponds more with the noise 
due to the subject's motion, since it has a low amplitude portion, which is followed by a high 
amplitude portion. Fig. 4 also shows that the second principal component 420 corresponds 

10 less with the noise due to the subject movement, since it does not show a distinct high 

amplitude portion. 

[37] Furthermore, a comparison of Figs. 6 and 7 shows that while the 
original mixture signals are not decorrelated, the principal components are decorrelated. Fig. 
6 shows a graph of the photocurrent at 655 nm vs. the photocurrent at 595 nm. As can be 
15 seen from this figure (Fig. 6), there is a wealth of mutual information between the two 
photocurrents, since the data from the two photocurrents appear to be aligned along the 
positively sloping diagonal line A-A. Fig. 7 shows a plot of principal component 2 vs. 
principal component 1. A review of this figure (Fig. 7) shows that the principal components 
are decorrelated, since there is no significant linear fit to the data in these coordinates. 
20 [38] In an alternate embodiment, before the processing according to step 

220, the time derivatives of the signals are obtained. For pulse oximetry, the third-order 
correlations of the pulse signals are often enhanced by taking the time derivative of the 
signals before performing PCA/ICA analyses. 

[39] Having decorrelated the data, the principal components are further 
25 processed by ICA processing to determine the independent components describing the 

photocurrent data (step 230). As described above, the principal components are processed to 
obtain a plurality of independent components, wherein a matrix of the plurality of signals 
corresponds with a matrix product of a matrix of the plurality of independent components and 
a matrix of mixing coefficients. In one embodiment, the independent components are 
30 determined by decorrelating the data by maximizing the sum of squares of the data set's 

higher-order cumulants. Since the skew of the time-derivative of the pulse signal is generally 
much greater in magnitude than that of interference, the ICA performance is enhanced by 
using a "contrast" function that was derived from the third-order cumulants of the derivatives 
of the signals. As stated above, a contrast function is any non-linear function which is 
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invariant to permutation and scaling matrices, and attains its minimum value in 
correspondence of the mutual independence among the output components. In an alternate 
embodiment, the independent components are determined by minimizing the estimated 
mutual information. Fig. 5 shows a plot 500 of the two independent components. As can be 
seen from this figure (Fig. 5), a first independent component 510 corresponds more with the 
noise due to the subject's motion, since it has a low amplitude portion, which is followed by a 
high amplitude portion. Fig. 5 also shows that a second principal component 520 apparently 
corresponds to a pulse component. Furthermore, Fig. 8 shows a plot of independent 
component 2 vs. independent component 1 . As can be seen from this figure (Fig. 8), the plot 
of independent component 2 vs. independent component 1 lies along the horizontal line at 
(independent component = 0) B-B, showing that the data sets have a minimal amount of 
mutual information, and thus can be approximated as independent data sets. 

[40] The decomposed data set of two independent components is further 
processed as follows. In one embodiment, sat is obtained from a ratio of mixing coefficients 
(step 240). In an alternate embodiment, one of the independent components is further 
processed to obtain the plethysmography signal (step 250), and the other independent signal 
is recognized as a measure of the interference signal (step 260). 

[41] Alternately, instead of, or in addition to taking signals that are 
measured at different wavelengths (step 210), signals are also obtained that are additions of 
signals from different times, thus an alternate embodiment of the present invention starts by 
processing signals obtained from mixing signals in time. 

[42] In pulse oximetry, embodiments of the present invention have the 
advantage of extracting the plethysmography signal in the presence of large motion 
interference and especially without requiring prior knowledge of saturation or pulse rate. 
Additionally, the present method is extendible to measurement of other physiological 
variables, such as pulse rate, blood pressure, temperature, or any other physiological variable 
measurable by non-invasive measurement of photocurrent provided an optical-based sensor. 

[43] Accordingly, as will be understood by those of skill in the art, the 
present invention which is related to blind source separation of pulse oximetry signals, may 
be embodied in other specific forms without departing from the essential characteristics 
thereof. Accordingly, the foregoing disclosure is intended to be illustrative, but not limiting, 
of the scope of the invention, which is set forth in the following claims. 
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